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ABSTRACT 

herween the shields, and conduction through the gas. One model is ba 

with 30, 45, 60, and 75 ML1 layers 

INTRODUCTION 

-ess s “2 itsxtts- 

liquids or to raise the system temperature. 

m a u ua aluminum or gold The mu, ..layer insulation systems - be used n bo* 
high-temperature and eugenic applications, 

arrangement, and heat transfer character.st.es are qu,.e d.fferent A detaded rev,ew 
MLI is provided by Tein and Cunnington’ 
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the spacer material becomes an issue at low imized mu systems for this type 

inner MLI layers isulating a cryogenic u greater distance than the warm outer 

of application, colder inner layers are separated by g referred as 

layers where radiation dominates the heat transfer. This type o mu 
V ariable Density because the layer spacing varies across the MLI system. The 3 

“1 system is controlled by the addition of bumper ^*““5 
Bumper strips are constructed by folding the dacron netting. Thus, thickness easily 

varied by varying the number of folds in the bumper strips 

111 

cryogenic applications. The M ( Y included several unique features 

XLn""roll wrap" inflation process provider l a robust “ 
assembly man-hours and seam heat leak A detailed description of the MHTB and 
MLI was provided by Hastings and Martin. 

The orbit hold test ,o“, 

: • «* -- “ 

"d ;lr~, f a consLriZsiX. ("“td p^^using the 

do however, tend to converge at the lowest boundary temperature of 160 K due to the 
decreased significance of radiation exchange between layers. 

lumpt" 
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Therefore, the focus in this paper is on analyttcal modeling of the MLI performance dunng 
orbital coast periods 

modeling overview 

orbital coast periods is discussed in subsequent sections. 

Layer-by-Layer Model 

The Laver-by-Layer model is based on McIntosh's work on a separated mode equation to 
materials. The total heat flux through the MLI is given by 


C|total Cjradialion 


+ 


q g as conduction + Toltd conduction 


( 1 ) 


The radiation heat transfer is, 


- o(T. 4 -Tc < )/(lte- +| fe- 1 ) 


( 2 ) 


where 


o = Stephan-Boltzmann constant = 5.675 E-8 W/m 2 -K 4 
T w = temperature of warm surface, K 
Tc = temperature of the cold surface, K 

e and 8c are the emissivities of the warm and cold surfaces, respec ive y. 


The gas conduction equation is provided by McIntosh who adopted an earlier formulation 
by Corruccini 6 


Qgas conduction 


= Ci P a (T w - T c ) 


( 3 ) 


where 

k g = Ci Pet 2 

kg = gas conduction, W/m -K 
P = gas pressure, pa 

C, = [(y + 1 V(t — 1)1 [R/8ttMT] ,,: 

a = accommodation coefficient 

y = Cp/Cv 

R = gas constant, 8 314 kJ/mol-K 
M = molecular weight of gas, kg/mol 
T = temperature of vacuum gage, normally 300 K 
For air, C, = 1 1666 and for helium 2.0998. 


Conduction through the solid is expressed as. 
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(4) 


C[so1id conduction 


= K, (Th-Tc) 


where. 


K s = C 2 f k/DX 
C 2 = an empirical constant 

f = relative density of the separator compared to solid material 

k = separator material conductivity, W/m-K 

DX = actual thickness separator between reflectors, m 

Curve fit equations have 

and glass paper as a function of temperature, 
provided: 


k = 0.01 7 + 7E-6*(800 - T) + 0.0228 ln(T) 
C 2 = 0.008 


(5) 


Modified Lockheed Model 

The Lockheed mode, also allows three 

conductivity of the spacer fibers is described as 

„ = ArN*') n T m (T h - T c )/N 

where A is an empirical coefficient, N* represents layer density £££££ 
average temperature S-J "represents the 

^ii^ tT£ conductive hea, transfer under free-molecule 

conditions, is shown as 

q„— = P (Y+ W- l)/(R/8nMT m )0.5P(T„-T c ) (7) 

perforated shields is given by 

q™inu«i = B e o (Th 4 67 " T C 4 ' 67 )/(N, + 1) (8) 

s is emissivity of the shields and B is an empirical parameter The total heat transfer can 
then be expressed as 

(9) 


C]ioul = C|solid conduction + Qradtation + Qgas conduction 
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Qtoui 


= A (N*) n T m (T h - Tc)/(N. +1) + B a (TV 


4.67 _ Tc 4 -6 7 )/(N 3 +1) 


+ CP(xJ)(TH (mM, -Tc^ ,) /(N s + D (1°) 

where P(x T) is the pressure within the insulation as a function of position and local 
teinperaturZ and coefficient A, B, and C as well as ,he exponents m and n are denved 
from the particular insulation system and intestial gas. 

For perforated aluminized shields, nitrogen gas, and glasstissue spacer material, the 
suggested Lockheed equation for the total heat transfer becomes the following. 

- 7 v 10' 8 fN*f 265 T (Th-T c )/N s + 7.07 x 10' 10 8 (Th 467 - Tc )/N« 
qiou.-7.30x 10 (N T.(l„ tc) + 1.46 x 10 4 P (Th 052 - Tc 0!2 )/N, (11) 

w ” ere q IOU i = total heat transfer in w/m 2 

s = 0.043 

N s = number of shields 


P = pressure in torr 
T = temperature in degree K 
N* = layers/cm 


In the original Lockheed equation, the spacer material and the shield hole sizes were 
different than those of the test article, therefore, the parameters A and B were adjusted fo 
the MHTB insulation The coefficient B in the Lockheed ^eqwuion accounts ;for ffie 
radiation between the shields and was provided for perforated shields of 0 ( . 

inches) hole diameter and a fractional open area of 0.0, To adj 

perforated shields used in the MHTB MLi [hole diameter of 1.27 cm (0^5 inches) and 
fractional open area of 0.02], reference 7 was used. The adjusted value of B is calculated 
fo be 4 944E-10 On the other hand, the coefficient A influences the conduction through 
the spacer The spacer material in the original Lockheed equation was tissuglass while in 
ffie es. article dacron spacer material was used. To modify the solid conduction term 
the dacron conductivity (unction provided by McIntosh is mcorporMK I into the 
conduction term. Therefore, the Modified Lockheed equation becomes the following. 


Q total = 2.4x10 


i -4 


(0.017 + 7E-6(800 - T) + 0.0228ln(T))* (N*) 

' - - - — -A — 0.52 


4.944 x 10’ 10 e (Th 467 


T c 467 )/N,+ 1.46x 10 4 P(T h ° -Tc 


2 63 (T h - Tc)/N, + 

0.52\ /xr (12) 


2 )/N s 


RESULTS AND DISCUSSIONS 


Figure 1 represents a schematic of the cryogenic tank with a foam/MLI combmation i with 
five different segments. The first segment is the foam. The second, third and fourth 
segments represent the three segments of MLI with different layer densit y and 
shields The layer density N*, and number of shields for MLI segments ru , , 

and 16 layers/cm and 10, 15, 20 shields respectively. Finally, the last legmen 
considered to be the shroud with an emissivity of 0.04. For the steady •state conditions, the 
heat flux through the MLI system is modeled using the modified Lockheed equation. The 
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model was comprised of described five different segments. Therefore, the heat transfer 
mode consisted of conduction through the foam, 3 segments of MLI using ‘he Lockheed 
ruation for each segment, and radiation exchange between the shroud and the last ML 
shield The Cold Temperature Boundary condition at the interior of he foam ' “d 
temperature, respectively. Experimental data «-nt^ned la Grid Temperature ^ 

the Hot Boundary Temperatures are 164 K, 235 K, and 
conditions, the following relations can be written as 

Qfoam = Qlsyerl = Qlayer2 = Qlayer3 = Qshroud ^ ^ 

TheO’s T, T 2 and T 3 are unknowns. Using equations ( 13 ), a system of four equations 
I,th four unknowns can be developed. Because of the radiation terms, this system of 
equations is nonlinear and is solved iteratively. First, for prescribed co d and ho. boundary 
temperatures initial MLI segment interface temperatures are assigned^ Then, heat rate 
through each segment is calculated. Using the heat rate and guessed tem P erat ^ e ’ 
resistant of each segment is computed. Then, new temperatures are computed using the 
nevv ^eat^atT and ^e si stance of each segment. The process is repeated until the so.u bon 
converges The convergence of the solutions implies that the total heat transfer rat 
through each segment is the same and the temperature of each MLI segment interface 
two consecutive iteration steps becomes equal within an allowable diflfere 
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Figure 1 Schematic of Foam/MLI System 

In the Layer by Layer approach, the foam, each shield layer, and the shroud are considered 
to be separate nodes. Therefore, the foam/MLI system model is comprised of 47 note 
For the steady state, heat flux through the MLI system was calculated using McIntosh 
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model and applying similar computational procedure as that of Lockheed model as 
described earlier. Based on the described methods, two spreadsheet modes w 
developed These spreadsheets could be useful tools in design optimiza ion a 
performance evaluation of MLI systems. 

Comparison between the heat flux values predicted by both models and those of the 

measured data is shown in Figure 2. Compared with the MHTB test data, for higher value 

of Hot Boundary Temperature, the heat flux values predicted by 

modified Lockheed equation are within 5% and 8%, respectively. Wh 

values of Hot Boundary Temperature, the predictions are within 34/o and 30/o for the 

Layer-by-Layer and Modified Lockheed, respectively. 

Comparison Between two models are extended in three hypothetical foam/MLI systems 
with die different shield layer densities than that of the test article. The MLI systems 
simulated are 30. 60, and 75 layers Each of these foam/MLI system is considered to hav 
similar physical and thermal properties as those of 45 layer foam/MLI system Figure 
4 and 5 illustarte the predicted heat flux for the MLI systems of 30. 60 and 75 layers 
respectively The predicted heat flux values at the highest Hot Boundary Temperature are 
within 6% As the Hot Boundary Temperature is lowered, heat fluxes predicted by o 
models converge to similar values The results for the foam/MLI system indicate that bo h 
the Layer-by-Layer and Modified Lockheed models simulations are consistently within 

6 %. 

In conclusion, either spreadsheet model can be utilized to predict the performance of 
either the MHTB type foam/MLI combination or the MLI alone. In the developmen o 
the analytical model for the MLI, it was observed that the larger, more widely spaced 
ventilation holes provided a significant radiation blockage advantage as compared with the 
standard closely space holes The variable density contributed to a weight reduction or 
performance increase due to reduced conduction with fewer layers to perform the same 
task The lack of seams, butt joints and structural support pins no doubt contributed to the 
performance improvement as well. The foam insulation not only enabled the elimination of 
a helium purge system, but also reduced the ground hold heat leak sufficiently to improve 
the effective density of the loaded hydrogen as compared with a “MLI” only concept. 
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Hot Boundary Temperature (K) 


Figure 2. Heat Flux Predictions for 45 Layers MLI 



Hot Boundary Temperature (K) 


Figure 3 Heat Flux Predictions for 30 Layers MLI 




Hot Boundary Temperature ( Degree K) 


Figure 4. Heat Flux Predictions for 60 Layers MLI 
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Figure 5 Heat Flux Predictions for 75 Layers MLI 






ACRONYMS AND ABBREVIATIONS 


cm 

centimeter 

hr 

Hour 

Kn 

Knudsen number 

K 

Degree Kelvin 

m 

Meter 

MHTB 

Multi-Purpose Hydrogen Test Bed 

MLI 

Multi-Layer Insulation 

VD-MLI 

Variable Density Multi-Layer Insulation 

W 

Watts 
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